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Abstract 

We construct new regular solutions in Einstein- Yang- Mills theory. 
They are static, axially symmetric and asymptotically flat. They are 
characterized by a pair of integers (k,n), where k is related to the 
polar angle and n to the azimuthal angle. The known spherically and 
axially symmetric EYM solutions have k = 1. For k > 1 new solutions 
arise, which form two branches. They exist above a minimal value of 
n, that increases with k. The solutions on the lower mass branch are 
related to certain solutions of Einstein- Yang-Mills-Higgs theory, where 
the nodes of the Higgs field form rings. 
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1 Introduction 



The well-known regular Bartnik-McKinnon (BM) solutions pQ and the corre- 
sponding non-Abelian black hole solutions ;2J, are asymptotically flat, static 
spherically symmetric solutions of SU(2) Einstein- Yang-Mills (EYM) the- 
ory They are unstable solutions, sphalerons t 3j, and are characterized by 
the number of nodes of the gauge field. Besides the BM solutions there are 
also asymptotically flat, static, only axially symmetric regular and black hole 
solutions jlj. These are characterized by two integers, the node number of 
their gauge field function(s), and the winding number with respect to the az- 
imuthal angle, denoted n. The spherically symmetric solutions have winding 
number n — 1, while the axially symmetric solutions have winding number 
n > 1. 

In SU(2) Einstein- Yang-Mills-Higgs (EYMH) theory, with a triplet Higgs 
field, gravitating monopole solutions arise t 5J. While gravitating monopoles 
with unit magnetic charge are spherically symmetric, the known gravitating 
multimonopoles possess axial symmetry [UJ. As in EYM theory, these EYMH 
solutions are characterized by two integers, the node number of the gauge field 
and the azimuthal winding number n, which corresponds to the topological 
charge of the monopoles. 

But EYMH theory allows for further static axially symmetric solutions, 
representing gravitating monopole-antimonopole pair, chain and vortex solu- 
tions [3 IHJ El- These solutions can be characterized by the azimuthal winding 
number n, and by a second integer m, related to the polar angle. For the 
monopole-antimonopole chains, which in flat space arise for n = 1 and 2, the 
integer m corresponds to the number of nodes of the Higgs field (and thus 
the number of poles on the symmetry axis) In vortex solutions, on the 
other hand, which in flat space arise for winding number n > 2, the Higgs 
field vanishes (for even m) on m/2 rings centered around the symmetry axis 

m 

The existence of monopole-antimonopole pair, chain and vortex solutions 
in EYMH theory |H1 E] immediately leads to the question, whether there 
might be analogous solutions in EYM theory, even though there would be no 
Higgs field participating in the subtle interplay of attraction and repulsion. 
In this letter we report the existence of one such new type of solution, related 
to vortex solutions in EYMH theory. 

In section II we present the EYM action, the axially symmetric ansatz 
and the boundary conditions. In section III we discuss the properties of the 
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new axially symmetric solutions, and we present our conclusions in section 
IV. 



2 Action and Ansatz 

We consider the SU(2) EYM action 

s = /(i4->(^)« (1) 

with Ricci scalar R, field strength tensor 

= d^Au - d v Ap + ie [A^ A v \ , (2) 

gauge potential A^ = r a A^/2, and gravitational and Yang-Mills coupling 
constants G and e, respectively. Variation of the action (JTJ) with respect to 
the metric leads to the Einstein equations, variation with respect to the 
gauge potential A^ to the gauge field equations. 

In isotropic coordinates the static axially symmetric metric reads [3] 

ds> = -fdfi + ^ + ^ + l ^fl V , (3) 

where the metric functions /, m and / are functions of the coordinates r and 
9, only. The z-axis (9 = 0, 7r) represents the symmetry axis. Regularity on 
the z-axis requires m = I there. 

For the gauge field we employ the ansatz [U |HJ |U] 



A-.dx^ 1 — 

M 2er 



r" {H x dr + (1 - H 2 ) rd9) - n (r r r l > k H 3 + t^HA r sin 



(1) 

Here the symbols r"' fc , Tq' and r™ denote the dot products of the cartesian 
vector of Pauli matrices, f — (r x ,T y ,T z ), with the spatial unit vectors 

e r n,h = (sin k9 cosrup, sin k9smmp, cos k9) , 
e 9 n ' k = (cos k9 cos rup, cos k9 sinrup, — sin k9) , 
ej 1 = (— sinrup, cosrup, 0) , (5) 

respectively. The gauge field functions Hi, i — 1 — 4, depend on the coordi- 
nates r and 9, only. For k = n = 1 and Hi = H 3 = 0, H 2 = 1 — H^ = w(r) 
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the BM solutions PQ are recovered, while for k = 1, n > 1, one obtains the 
axially symmetric solutions of [3]. The new solutions reported here are ob- 
tained for k > 1. They are related to EYMH solutions with m = 2k in the 
limit of vanishing Higgs field . 

The ansatz is form-invariant under the abelian gauge transformation [3] 

U = exp (~r£r(r,0)) . (6) 

We fix the gauge by choosing the gauge condition |U |HJ E] 

rd r Hx - d e H 2 = . (7) 

To obtain asymptotically flat solutions which are globally regular and 
possess the proper symmetries, we need to impose appropriate boundary 
conditions [H |BJ E] ■ At the origin we impose the boundary conditions 

d r f = d r m = d r l = , Hx = H 3 = H A = , H 2 = 1 , (8) 

at infinity we impose 

/ = m = l = l , Hx = H 3 = , H 2 = 1 - 2k , H 4 = 2 am(k0) / sin 9 , (9) 

and on the z-axis we impose 

d e f = d e m = d e l = 0, Hx = H 3 = 0, d e H 2 = d e H 4 = . (10) 

We further introduce the dimensionless coordinate x, and the dimension- 
less mass /i, 

x = . r , u = — M . (11) 

3 Numerical Results 

Subject to the above boundary conditions, we solve the system of seven 
coupled non-linear partial differential equations numerically. To map spatial 
infinity to the finite value x — 1, we employ the radial coordinate jU |H1 E] 

^ = T^-- (12) 

1 + X 
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The numerical calculations are based on the Newton-Raphson method, and 
are performed with help of the program FIDISOL ^Uj. The equations are 
discretized on a non-equidistant grid in x and 6. Typical grids used have 
sizes 70 x 30, covering the integration region < x < 1 and < 9 < ir/2. 
For the method, it is essential to have a good first guess, to start the iteration 
procedure. For the k = 1 EYM solutions, the n = 1 BM solutions serve as a 
first guess, and then the 'parameter' n is varied (via unphysical noninteger 
values) to obtain solutions with higher winding number n \\}. For the k = 2 
solutions we employ the m = 4, n = 4 EYMH vortex solution with (almost) 
vanishing Higgs field as a first guess j^j, and then again vary n. Similarly, for 
the k = 3 solutions we start from the m = 6, n = 6 EYMH vortex solution. 

In EYMH theory, the monopole-antimonopole pair (MAP) solution is 
obtained, when m — 2k — 2, n — 1. Here a monopole and an antimonopole 
are located symmetrically on the z-axis. Similarly, one might try to find a 
sphaleron-antisphaleron pair (SAP) solution in EYM theory. The boundary 
conditions of such a SAP solution, however, do not differ from those of a 
BM solution. Consequently, when the EYMH coupling constant a = y/AnGv 
(where v is the Higgs field expectation value) is varied, the gravitating MAP 
solution starts from the flat space solution, reaches a critical solution at a 
maximal value of a, and then with decreasing a smoothly reaches the BM 
solution in the limit a — > (after rescaling of the radial coordinate) [Jj. 
Consequently, we do not find a SAP solution in EYM theory, and similarly, 
we do not find new solutions, when m = 2k = 2, n > 1, since their boundary 
conditions agree with those of the axially symmetric solutions of jl]. 

When m = 2k = 4, however, the boundary conditions eq. © for the 
gauge field differ from those of the known solutions Ej. We therefore 
may expect the existence of new EYM solutions, subject to these boundary 
conditions. Indeed, for m = 4, n = 4, for instance, when a is varied again, 
the gravitating EYMH vortex solution starts from the flat space solution [S] , 
reaches a critical solution at a maximal value of a, and then with decreasing 
a smoothly reaches a new EYM solution in the limit a — > (after rescaling 
of the radial coordinate) jU]. 

Let us in the following refer to the solutions characterized by the integers 
k and n as (k, n) solutions. In Fig. la we exhibit the energy density of the 
new (2,4) EYM solution. Its energy density has a pronounced maximum 
on the p-axis, thus the energy density is torus-like. The energy density of 
this (2, 4) solution is thus similar to the energy density of the (lowest mass) 
known (1,4) EYM solution exhibited in Fig. lb. But the maximum of 
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the (2,4) solution is smaller than the maximum of the (1,4) solution, and it 
is shifted slightly inwards. In addition, the (2, 4) solution has a small saddle 
point on the p-axis, located further outwards. On the other hand comparison 
with the corresponding EYMH solution [8 j shows, that the energy density 
of the EYM solution differs significantly from the energy density of the flat 
space (2, 4) EYMH solution, whose energy density exhibits two tori, located 
symmetrically with respect to the xy-pld,ne |8,. 

e k=2, n=4, lower branch 

0.12 r 




Fig. 1 The energy density of the (2, 4) EYM solution on the lower mass branch 
(a) and of the (1,4) EYM solution (b). 

The existence of a k = 2 EYM solution with winding number n = 4, 
immediately suggests that analogous solutions with winding number ti ^ 4 
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should exist as well. Indeed, by increasing n, a whole set of such new solutions 
is obtained numerically, one for each integer n. On the other hand, a decrease 
of n does not lead to new solutions. Instead a bifurcation at the non-integer 
and thus unphysical value of n — 3.7175 is obtained, where a second branch of 
solutions with higher mass appears. Consequently there is a second solution 
for n — 4, and also for each higher integer n, but there are no such solutions 
for integer n < 3. This is illustrated in Fig. 2a, where we exhibit the mass 
of these (2, n) EYM solutions. Only the integer values of n in the figure 
correspond to physical solutions. In Fig. 2b we illustrate the value of the 
metric function / at the origin. We note, that for (l,n) solutions (with a 
single node) we find only one branch, which is also exhibited in Figs. 2. 



(a) 




(b) 




Fig. 2 The mass (a) and the value of the metric function / at the origin (b) of 
the (1, n), (2, n) and (3, n) EYM solutions is shown as a function of n. 



To illustrate the dependence of the solutions on n, we exhibit the energy 
density of the (2, n) solutions with n = 4, 6, 8 on the lower (mass) branch in 
Fig. 3a and on the upper branch in Fig. 3b. With increasing n the maximum 
of the energy density decreases and moves outwards for the solutions on both 
branches. For the same n, the maximum of the upper mass solutions is higher 
and located further inwards than the maximum of the lower mass solutions. 
The small saddle point turns into a small maximum on the upper branch, 
thus the energy density of the solutions on the upper branch has a more 
complicated structure as seen in Fig. 4, where a surface of constant energy 
density is shown for the (2,4) solution. We exhibit the gauge field function 
Hi for the same set of solutions in Figs. 5. 
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k=2, lower branch 



k=2, upper branch 
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(b) 
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Fig. 3 The energy density of the (2, n) EYM solution on the lower mass branch 
(a) and on the upper mass branch (b) for n = 4, 6, 8. 



k=2,n=4, upper branch 



Si 




Fig. 4 A surface of constant energy density (e = 0.005) of the (2, 4) EYM solution 
on the upper mass branch. 



For k = 3 we obtain a similar pattern of solutions. But the new type 
of solutions appears only for n > 5. The mass and the value of the metric 
function / at the origin of these (3, n) solutions are also shown in Figs. 2. 
The energy density for (3, n) solutions with n = 6,8, 10 and the gauge field 
function Hi are shown in Figs. 6 and 7, respectively. Based on these results, 
we conjecture the existence of (k, n) solutions with k > 3, but expect that 
they will appear only for still larger values of n. 
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k-2, lower branch 



k-2, upper branch 
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(b) 



-1 -0.5 0.5 1 1.5 2 2.5 3 

iog 10 M 



Fig. 5 The gauge field function H2 of the (2, n) EYM solutions on the lower mass 
branch (a) and on the upper mass branch (b) for n = 4, 6, 8. 
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Fig. 6 The energy density of the (3,n) EYM solution on the lower mass branch 
(a) and on the upper mass branch (b) for n = 6, 8, 10. 
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Fig. 7 The gauge field function H2 of the (3, n) EYM solutions on the lower mass 
branch (a) and on the upper mass branch (b) for n = 6, 8, 10. 
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4 Conclusions 



We have constructed numerically new static regular solutions of EYM theory 
with axial symmetry. These solutions are characterized by the integers (k, n), 
related to the azimuthal and polar angles, respectively. In particular, we have 
obtained solutions for k = 2, n = 4 — 10, and k = 3, n = 6 — 10. 

Like the previously known (l,n) solutions [4J, the (2,n) and (3, n) solu- 
tions on the lower branch have a torus-like shape. On the upper branch of 
the (3, n) solutions double tori appear for larger n. The (1, n) solutions (with 
one node) most likely exist for any integer n > 1 and form a single branch. 
The new (2,n) and (3, n) solutions have lower bounds on n, n = 2k, which 
imply the existence of two branches of solutions for both values of k. We ex- 
pect, that the (2,n) and (3,n) solutions represent only the first sequences of 
new solutions, and conjecture the existence of (k, n) solutions also for higher 
values of k. 

The k = 1 BM solutions are unstable |3|, and there is all reason to be- 
lieve, that the new (k, n) solutions are also unstable. Having constructed 
the regular axially symmetric solutions it appears straightforward to con- 
struct analogous black hole solutions, and, in particular, to look for rotating 
solutions [T2] . 
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